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Предлагается новый алгоритм идентификации параметров эквивалентной схемы замещения электрического разряда. Подход 
основан на решении интегрального уравнения I рода относительно функции переходной проводимости, по которой затем про- 
исходит определение параметров схемы замещения. Использование сглаживающих сплайнов и оригинального регуляризирую- 
щего алгоритма, учитывающего погрешности задания ядра интегрального уравнения позволило получить устойчивый алгоритм 
идентификации параметров. Проведенные исследования алгоритма показали высокую вычислительную эффективность и хоро- 
шую точность идентификации параметров. 


1. Введение 

Одно из наиболее интересных с физической 
точки зрения и практически важных направлений в 
различных областях техники является барьерный 
разряд. В частности барьерный разряд использует- 
ся для очистки воды, плазменных технологий, 
травления и т. д. Однако сильная пространственная 
неоднородность и малая длительность физических 
процессов, протекающих в барьерном разряде, 
крайне затрудняет изучение этого явления. 

При феноменологическом описании электри- 
ческих разрядов (в частности и барьерного разряда) 
широко используется их описание как объектов 
электрической цепи [1]. В основе таких подходов 
лежит замена электрофизических явлений, проис- 
ходящих в разряде, явлениями, происходящими в 
электрической цепи, состоящей из определенных 
электрических элементов (сопротивления, емко- 
сти, индуктивности). Такую электрическую цепь 
будем называть эквивалентной схемой замещения 
электрического разряда. 

При исследовании физики разрядов доступны- 
ми для измерения являются напряжение 11(1) и ток 
1(1) в цепи с разрядным промежутком. Поэтому 
возникает задача определения параметров эквива- 
лентной схемы замещения по зарегистрированным 
значениям функций 11(1), І({). По сути, имеем зада- 
чу идентификации параметров эквивалентной схе- 
мы замещения. 

В работах [2, 3] значения параметров находятся 
по переходной проводимости §(1). В свою очередь 
функция §(1) определяется как решение интеграль- 
ного уравнения-свертки Вольтера I рода, что явля- 
ется некорректно поставленной задачей [4, 5]. Од- 
нако в этих работах не учитывается ряд важных мо- 
ментов, связанных с решением этой некорректно 
поставленной задачи, что отрицательно сказалось 
на точности идентификации параметров. К таким 
моментам можно отнести: недостаточная устойчи- 
вость используемого алгоритма дифференцируемо- 
го зашумленной функции Ѵ(і) при вычислении 
ядра интегрального уравнения; не учет случайной 
погрешности задания ядра уравнения, как на этапе 


построения регуляризированного решения, так и 
при выборе параметра регуляризации. 

Поэтому в данной работе предполагается устой- 
чивый алгоритм идентификации параметров схемы 
замещения, основанный на регуляризирующем ал- 
горитме решения интегрального уравнения-сверт- 
ки с неточно заданным ядром [6] и в полной мере 
учитывающей вышеназванные моменты. 


2. Постановка задачи 


Если действующее в цепи напряжение имеет 
импульсную форму, то переходный процесс, про- 
исходящий в разрядном промежутке в терминах 
11(1) и 1(1) удобно описать с помощью интеграла 
Дюамеля [7]. Если действующее напряжение 11(1) 
является финитной функцией, т. е. вне интервала 
[О, Гу] обращается в нуль, то для тока в цепи спра- 
ведливо выражение 


Щ = Щ0)ш«) + 


I 


сііі(т) 

СІТ 


§(1-т)Фт, 


где §(1) - переходная проводимость. Как правило, 
значение Щ 0)=0 и поэтому приходим к интеграль- 
ному уравнению-свертки Вольтерра I рода: 


ДО = _[&(* -г) 

о 


сІт 


(1) 


Функцию §(1—т) называют ядром интегрального 
уравнения, 1(1) - правой частью уравнения. 

Интегральное уравнение (1) необходимо ре- 
шить относительно функции §(1), что является не- 
корректно поставленной задачей, а затем по функ- 
ции §(1) определить параметры эквивалентной схе- 
мы замещения электрического разряда. 

Таким образом, задача идентификации параме- 
тров эквивалентной схемы замещения включает 
следующие этапы: 

Этап 1 . Вычисление производной по изме- 

сіт 

ренным (с погрешностями) значениям функции 11(1). 

Этап 2. Решение интегрального уравнения (1) 
относительно функции §(1). 
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Этап 3. Определение (возможно и по виду 
функции §(!)) структуры эквивалентной схемы за- 
мещения и параметризации функции §(1). 

Этап 4. Оценивание параметров функции §(1) и 
вычисление по этим оценкам величин сопротивле- 
ний, емкостей и индуктивностей, входя щ их в экви- 
валентную схему замещения. 

Решение сформулированной задачи идентифи- 
кации будем рассматривать при следующих пред- 
положениях: 

1. Функция 11(1) отлична от нуля на интервале 
(О, Гу] (т. е. является финитной) и измеряется на 
этом интервале в моменты /]=А-у,у'= 0,1,...,Л^-1, где 
7Уу=еЩ[Гу/Д]+1, А - шаг дискретизации, еп([г] - 
целая часть вещественного числа і. Измеренные 
значения Ц допускают представление 

С , - 17( /А) і 4’/. ] = 0,1,..., Ы и -1, 

где 4) - случайные величины с математическим 
ожиданием М(^)= 0, дисперсией Д4/)=с>/ и ото- 
бражающие погрешности измерения напряжения. 

2. Функция §(1) отличается от нуля на интервале 
[0,Г 8 ]. 

При этих предположениях функция 1(1) являет- 
ся финитной с интервалом определения [0,7}], где 
Ті=Т п +Т е . 

3. Функция 1(1) измеряется на интервале [0, 7}] в 
моменты (=Ау, у— 0,1,..., Л}— 1, где Л}=епІ[7}/Д]+1. 
Измеренные значения допускают представление 

/>/(/,) + ц ; , 7 = 0,1,...,^ -1, 

где ц 7 - случайные величины с числовыми характе- 
ристиками Л/(77 7 )=0, Д77 7 )=(5, 2 . 

Кратко остановимся на алгоритмах решения 
каждого этапа рассматриваемой задачи идентифика- 
ции параметров эквивалентной схемы замещения. 

3. Вычисление производной 

по измеренным значениям напряжения 

Ядром интегрального уравнения (2) является 
сШ(і) ти ч тт 

производная напряжения і!(1). Известно, 

Ат 

что операция дифференцирования является некор- 
ректно поставленной задачей (в частности малые 
ошибки могут вызвать сколь угодно большие 
ошибки в производной). 

Для устойчивого дифференцирования функции 
Щі), заданной измеренными в моменты значения- 
ми Ѵ(1) в качестве приближения для 11(1) примем 
сглаживающий кубический сплайн (СКС) А я (/). На- 
помним, что кубическим сглаживающим сплайном 
называется кубический полином, удовлетворяю- 
щий условиям: 

1. На каждом интервале [7 7 ,( +1 ] 8 Х (1) имеет сле- 
дующее представление 

(0 = а ] + Ь / (I - Г. ) + с ] (1-1Щ 2 + А ) (I - 1; )\ 
где (<!<(„. 


2. Функция А я (7) имеет непрерывную вторую 
производную на всем отрезке [0,7},]. 

Вычисление коэффициентов о,, Ь р с,, ф СКС 
(которые зависят от параметров сглаживания Я) по- 
дробно изложено в работах [5, 8] и для их однознач- 
ного вычисления примем краевые условия вида: 

Д'(0) = 0; Д'(Г с/ )=0. (2) 

Эти условия соответствуют типичной форме 
импульса напряжения 11(1) (см. рисунок). Можно 
показать, что СКС с условиями (2) доставляет ми- 
нимум функционалу 

Т и N„-1 

[(Г(0) 2 л+я-Х^.(/д.)-с.) 2 

о 7=0 

среди всех функций Д/) с интегрируемым квадра- 
том второй производной и удовлетворяющих усло- 
вию (2). 

После вычисления коэффициентов СКС пер- 
вую производную 8[(і) (являющуюся оценкой для 

производной ) МОЖНО вычислить по формуле 
Ат 

я/ (0 = Ъ і + 2 С і (I - ?,. ) + 3 А } (1-1^-, (3) 

где (< /</}+,. 

Основной сложностью при построении СКС яв- 
ляется выбор параметра сглаживания Я, который 
может изменяться в пределах от 0 (сглаживающий 
сплайн становиться интерполяционным, проходя- 
щим через значения Ц, т. е. >%(()= Ц) до оо (СКС ста- 
новиться прямой линией). Если Я окажется малым, 
то в сплайне будут присутствовать высокочастот- 
ные составляющие, обусловленные погрешностями 
4), которые будут особенно проявляться в производ- 
ных сплайна в виде высокочастотных осцилляций. 
Если этот параметр будут слишком большим, по- 
строенный сплайн окажется «переглаженный» и в 
нем будут сильно сглажены передний и задний 
фронты импульса 11(1), что отрицательно скажется 
на точности вычисления первой производной. 

Можно выделить два подхода к выбору параме- 
тра Я: оценивания Я из условия минимума среднек- 
вадратичной ошибкой сглаживания [5] и выбор Я 
по точностным характеристикам сплайна [9]. Оста- 
новимся на втором подходе, более подходящем для 
решаемой задачи дифференцирования 11(1). 

В этом подходе сглаживающий сплайн интер- 
претируется как выходной сигнал некоторого 
фильтра (сплайн-фильтра), на вход которого по- 
ступает дискретная последовательность, состоя- 
щая из измеренных значений 7} функции Щі). При 
такой трактовке сглаживающие свойства сплайна 
можно определить через его аппаратную функцию 
И х (1), которая характеризует систематическую 
ошибку сглаживания и дифференцирования: чем 
меньше «ширина» функции И х (1), тем меньше си- 
стематическая ошибка. В качестве числовой харак- 
теристики аппаратной функции принимается ее 
ширина А й (Я): 
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Рисунок. Вычисление производных функции і 1(1) 
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Физическая трактовка этой характеристики для 
задачи дифференцирования достаточно проста: в 
сглаживающем сплайне и его производной сохра- 
няются (с небольшими амплитудными искажения- 
ми) составляющие функции 11(1) и производной 
Щі), если их ширина больше ширины аппаратной 
функции И,(і). Задавая «предельный» размер А„ р со- 
ставляющих, которые должны сохраниться в 
сплайне, значение А можно определить из решения 
нелинейного уравнения: 

А*(Я) = А чі . (4) 


Аппаратная функция к,(і) вычисляется по фор- 
муле 





Н х (б))е""'сІ со. 


Частотная характеристика Ні со) сплайна, опре- 
деляется выражением [9] 

Н х (со) = 


2 

( ■ с°> Т ^ 
»»(— I 

2 

1 - соз(й)А) 

А со 2 

соТ 

V 2 


до + 2</, С 08 (й)А)+ 2 ді со$(2й)Д) 


2А 6 Хр А 4 А о Хр 
где до = — + — -ср = -- — - д2 = —-р.= р 

- весовые множители функционала, А - шаг дис- 
кретизации. 

Проиллюстрируем изложенный подход к выбору 
параметра сглаживания Я результатами следующего 
вычислительного эксперимента. На рисунке, а , по- 


казан график функции 100 11(1) (кривая 1) и ее «точ- 
ной» производной —11(1) (кривая 2). Значения 
сіі 

функции умножаются на 100 для того, чтобы в мас- 
штабе рисунка эта функция отличалась от нуля. 
Кривая 3 соответствует производной интегрального 
сплайна А 0 (?), построенного по измеренным (с по- 
грешностью) значениям Ц,і= 1,2,... ,#[,=240. Относи- 
тельный уровень погрешностей задавался равным 
0,05. На рисунке, б, показан график «точной» про- 
„ с іЬ (і) 

изводнои — — и значения производной, вычи- 
сіт 

сленной по интерполяционному кубическому 
сплайну. Видны значительные осцилляции этой про- 
изводной, характерные для дифференцирования не- 
точно заданных функций. На рисунке, б, приведен 
„ (Ш (I) 

график производной и значения производ- 

ит 

ной >% '(і), вычисленной по сглаживающему кубиче- 
скому сплайну 8 } (і), см. (3). Видно достаточно хоро- 
шее (по сравнению с производной интерполяцион- 
ного сплайна) совпадение этих производных. 


Параметр сглаживания выбирается из решения 
ур. (4) при А л/ ,=5- 10 3 с (интервал дискретизации 
Д, ф =2,5-10~ 4 с). Величина А„ р задавалась равной поло- 
вине ширины переднего фронта импульса напряже- 
ния Щі), что позволило в производной 8 } '(і) сохра- 
нить «тонкие» детали производной 11(1) (в частности 
на интервале [0, 0.01]). Небольшие колебания 81(1) 
производной на интервале [0.01, 0.03] будут на вто- 
ром этапе интерпретирования как погрешности за- 
дания ядра интегрального уравнения и будут учтены 
при построении регуляризированного решения это- 
го интегрального уравнения. Для этого представим 
значения производной сплайна 8 '(1) в узлах ( в виде: 


= 0 , 1 ,...,^ - 1 . 
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Случайные величины § отображают ошибки в 

вычислении производной по сглаживающему 

сплайну 8)(і). Если погрешности ^ измерения Щ) 
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